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ABSTRACT 



Context. We study the evolution and fate of solar composition super-massive stars in the mass range 60 — 1000 Mq. Our study 
is relevant for very massive objects observed in young stellar complexes as well as for super-massive stars that could potentially 
form through runaway stellar collisions. 

Aims. We predict the outcomes of stellar evolution by employing a mass-loss prescription that is consistent with the observed 
Hertzsprung-Russell Diagram location of the most massive stars. 

Methods. We compute a series of stellar models with an appropriately modified version of the Eggleton evolutionary code. 
Results. We find that super- massive stars with initial masses up to 1000 Mq end their lives as objects less massive than ~ 150 Mq . 
These objects are expected to collapse into black holes (with M ^ 70Mq) or explode as pair-instability supernovae. 
Conclusions. We argue that if ultraluminous X-ray sources (ULXs) contain intermediate-mass black holes, these are unlikely to 
be the result of runaway stellar collisions in the cores of young clusters. 
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1. Introduction 

In this paper, we study the structure and evolution of very 
massive stars (VMS), defined as objects with masses of 60 
up to 150Mo, as well as super-massive stars (SMS) with 
masses in the range 150 - 1000 Mq. 

The interest in the upper limit of stellar masses and 
the evolution and fate of the most massive stars in 



ultraluminous X-ray sources (ULX, see 


Fabbiano 19891 


Colbert & Milleil|2005llFabbiano & White 


I2OO6I; Fabbianol 


2006; ,Soria, 20061. and references therein). These objects 



are most commonly interpreted as binaries involving ei- 
ther sub-Eddington accretion onto an intermediate mass 
black hole (IMBH) with mass ~ (10^ - 10^) Mq or super- 
Eddington accretion onto a stellar mass black hole with 
mass ~ 10 Mq. In the latter case, beaming or support 
of super-Eddington luminosity by an accretion disk is re- 
quired. Currently, the issue of the black hole mas ses in 
ULXs remains unresolved (see, e.g., Fabbiano 20061 ). 

It has been argued that IMBHs may als o reside in the 
cores of some globular clust ers 



se e, e.g. 



Gerssen et aL. 

2OO2I I2OO3I: iGebhardt etHI l2005t IPatruno et al.l l2006t 



Feng fc Kaa"re^l2006l ). but see lBaumgardt et al.l (|2003allbl ) 
for counter-arguments. The masses of the putative black 
holes in the cores of well-known globular clusters such 
as 47 Tu c and NGC 6397, are expected to be ~ (10^ - 
10^) Mq (|De Riicke et al.ll2006l ). 

ULXs are observed in both spiral and elliptical galax- 
ies, i.e. in environments with diverse metal abundances 
and star formation rates. In this paper, we focus on 
stars with solar initial composition (X=0. 7,^=0. 02). 
This choice is partly motivated by the high metal- 
l icity of the s t arbur s t galaxy M82 - w it h Z :^ Zq 



( McLeod etakl Il993l : lOrigha et ai] 12004 iMavva et al" 



20061 ) - which contains one of the most promising can- 
didate ULXs, M82 X-1, and that may host a black hole 
of interr nediate mass. a.s arg ue d by e.g. iPtak fc Griffithsl 
(■I999I): iKaaret et all J2OO1I): iMatsumoto et al.l (I2OOII) 



Strohmaver fc Mushotzkv ( 20031 ). but see Okaiima et al 
( 20061 ) for arguments in favour of a stellar mass black hole 
in M82 X-1. 

The current observational est imate of the upper cut- 
off of stellar masses is ~ 1 50 Mq ("Masse y fc Hunter 1998; 
Weidner fc Kroupall2004t lFigeii l2005). However, for solar 
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chemical composition, it is expected that even such mas- 
sive stars rarely produce black holes more massive than 



2 



Yungelson et al.; Super-massive stars 



~ 20 Mq (|Maedej|l992l : iFrver k Kalogerall2001f) . due to 
copious mass loss during both the hydrogen and heUum 
burning stages. If ULXs and young globular clusters really 
harbour black holes with masses exceeding several 10 Mq, 
the problem of their formation becomes a challenge for the 
theory of stellar evolution with mass loss. 

Related intriguing problems involve the formation and 
evolution of the luminous stars found in the central par- 
sec of the Galaxy, e. g., the S-stars and the IRS 13 and 
IRS16 con glom erate s (jSchodel et allbOOSt ILu et al.|[2004 ; 
Maillard et al.ll200# l. or in the R136 stellar complex of 



the 30 Doradus nebula i n the Large Magellani c Cloud 
(|Walborn fc Blades! [iQQTt iMassev fc Huntei]ll998l ). 

Returning to the problem of IMBH formation, 
there are currently two models in favour. One in- 
volves the gradual accumulation of mass by accre- 
tion onto a seed black hole, which, while swallowing 
gas and stars in a stellar cluster, r nay grow to th e 



interm ediate mass range (see, e.g., iMiller fc Colbert 



2004 and references therein). The other scenario con- 
siders the collapse of an object which descends from 
an SMS formed by hierarchical runaway merger of or- 



or after core collapse (see, e.g.. 


— „ — , — 

Portegies Zwart et al. 


1999; .Portegies Zwart & McMillan 


2002: McMillan et al. 


2004; Giirkanetal. 2004; Portegies Zwart et al. 20041 


Freitag et all 2006). 



The motivation for the current study stems from the 
latter scenario. It assumes that stars are born with a 
wide distribution of masses ranging from the hydrogen- 
burning limit up to a maximum of about 150 Mq. More 
massive stars observed in the Galaxy may originate from 
stellar mergers. Most of these may be the result of coa- 
lescence of components of binaries in common envelopes. 
The mass of binary me r ger pr o ducts may be up to 3 00Mf7) . 



Portegies Zwart et al 



Portegies Zwart et al 



(1999): 



20041) 



Miller fc HamiltonI (|2002l) : 



Gtirkan et al.l (20041) have 



demonstrated by means of detailed N-body simulations, 
in which simple stellar evolution was taken into account, 
that there is a range of initial conditions where stellar co- 
agulation drives the mass of a star to ~ 1000 Mq. The 
resulting star burns-up quite quickly, and the process of 
hierarchical merging in the cluster core terminates as soon 
as the first massive stars experience supernovae and col- 
lapse into black holes. The collision runaway process ter- 
minates as the mass loss from the explosions of massive 
stars in the cluster center drives the expansion of the clus- 
ter core. This scenario was used to explain the large black 
hole mass of M82 X-1 (jPortegies Zwart et al.ll2004[ ). 

The studies of hierarchical mergers take into account 
the possibilities of coUisional mass loss, mass loss from 
stellar winds, and rejuvenation of merger products due to 



^ Note that a high stellar mass inferred from luminosity may 
be an artifact of unresolved binarity. For instance, IRS16 SW 
turned out to be a b inary with two almost identical compo- 
nents of ^ 50 Mq , (jDePov et all l2004l : iMartins et all l2006l : 
iPeeples et al]|2007t ). 



fresh hydrogen supply in collisions. However, stellar evo- 
lution is treated rather crudely in these simulations, using 
extrapolations by several orders of magnitude for stars 
that typically have zero-age main-sequence (ZAMS) mass 
< 100 A/q. The possible formation of core-halo configura- 
tions and the existence of an upper stellar mass limit are 
generally ignored. The treatment of mass loss in stellar 
winds is particularly uncertain. In principle, the winds of 
very massive and super-massive stars may be so strong 
that the cluster core expands in such a dramatic way that 
it stops being dominated by collisi ons, in which case the 



hierarchical me rger is terminated (jPortegies Zwart et al 
19991 : IVanbe^ren 2005) . 



We aim to refine evolutionary calculations for merger 
products, and as a first step, we study the evolution of 
solar composition VMS and SMS over the mass range 60- 
1000 Mq. We construct chemically homogeneous models 
for these stars and confirm the existence of their upper 
mass limit (~ 1000 Mq), and study their evolution with 
mass loss through the core hydrogen and core helium- 
burning stages. This allows us to determine the mass and 
nature of the pre-supernova objects and to predict the fate 
of these objects. 

The results of our evolutionary computations are pre- 
sented and discussed in Sect.^J a comparison with obser- 
vations is given in Sect. [3l the fate of SMS is discussed in 
Sect, m Conclusions of our study follow in Sect. [5l In the 
Appendix, we discuss a number of individual stars in close 
proximity to the Humphreys-Davidson (HD) limit. 

2. The stellar evolution calculations 

We compute the inner structure and evolution of single 
non-rotating SMS. Initially, all stars are chemically ho- 
mogeneous. The products of hierarchical merging are ex- 
pected to grow in mass in discrete steps by the injection 
of other stars and as a result the merger product may be 
rapidly rotating. In this paper, we ignore the effects as- 
sociated with merger events themselves, but concentrate 
on the ID (non-rotating) evolution of the final object, the 
super-massive star (SMS). 

Another open question, which has persisted for 
deca des, is that of the p ulsational stability of SMS (see, 
e.g., iBaraffe et al.ll200l[ ). Existing nonlinear calculations 
of the effect s of pulsation s on m ass loss differ by an order of 
magnitude. Appenzeller (1973) found that pulsationally- 
driven mass loss for solar metallicity stars with M > 
300 Mq occurs on a timescale that is shorter than the 
core-hydrogen burning timescale, but IPapaloizou ( 1973 ) 
claimed that the mass-loss rate is an order of magni- 
tude lower and evolutionarily insignificant. Wolf-Rayet 
(WR) stars and luminous blue variables (LBVs) are the 
descendants of main-sequence stars and are likely to be 
subje ct to pulsational instabilities (jPadevev fc Novikova 
20041 ) . Awaiting the resolution of this problem, we do not 
explicitly consider the possibility of vibrational mass loss 
in our calculations. 
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The evolutionary computations were carried out by 
means o f an appropriately modified version of the 
Eggleton (|l97ll . priv. comm. 2003) evolutio nary code. 
The input physics has been described in iPols et al. 
( 1995f ). but here we briefly list the sources of opacities 
and nuclear rates. The opacity tables ar e constructed 
using OPAL tables of llglesias &: Rogers! dlOOS '). low- 
tempe ratures opacities are from [Alexander fc Ferguson 



1994 • whilst iHubbard fc Lampd ITgegl ) and lltoh et al 



19831 ) list conductivities for degenerate matter. The nu- 
clear evolution of ^H, "iJe, ^^C, ^^Ne^^^Oj^J^e^^lMaj^ 



follow ed us ing the reaction rates fr om lCaughlan fc Fowler 
( 1988[ ) and Caughlan et al. ( 1985 ). The equation of state 
is based on the principle of Helmholtz free energy min- 
imisation and provides the physical quantities such as the 
pressure, density, specific heats etc. as a function of a pa- 
rameter related to the electron degeneracy and tempera- 
ture. 

An initial solar chemical composition (^=0.7 and 
Z—Q.Q2) is assumed, and the computations may be rel- 
evant to the Galactic objects such as the most massive 
members of Arches, Quintuplet . IRS13, IRS16, NGC3 603, 
Westerlund clusters (see, e.g.. IPortegies ZwartI 12004 for 
further references) and, in case they do exist, to SMS 
in those external galaxies that have chemical abundances 
comparable to the composition of the Milky Way. We con- 
sider stars with surface helium abundance Yc > 0.4 and 
Teff > 10000 K to be WR stars. 

2.1. Homogeneous supermassive stars 

With the opacity increasing outwards, a star develops 
a structure where most of the mass is concentrated in 
the compact convective core, whilst only a small frac- 
tion of the mass is located in a very extended radiativ e 
envelope ("core-halo configuration", iKatol 1X9851 Il98d ). 
More massive stars tend to be more extended, and the 
expansion of the envelope is caused by the increase of 
the Eddington luminosity Lsdd = ^TTcGMr/>c, where 
X is the flux-mean opacity in the radiative outer lay- 
ers of the star. The maximum mass for a stable star is 
reached when the luminosity of the star reaches L^dd 
in the photosphere. Using Compton scattering opacity 
X = 0.2 (1 + X)/{l + 2x IQ-^T), Kato succeeded in con- 
structing homogeneous hydrogen-rich solar metallicity 
stellar models of 10^ M0. The inclusion of the Kramers 
term in the opacity resulted in a reduction of the limit- 
ing Z AMS mass to w 10^ Mq (jMenshchikov fc Tutiikov 
These last authors also found an upper limit of 



ab out 2500 M cT, for h ehum-rich stars (X=0, Y=0.97). 
Ishii et al. ( 19991 ) explored the core- halo effect using 



modern OPAL opacities and found that, in this case, the 
upper mass limit for hydrogen-rich stars drops drastically. 
For stars of solar composition, this limit is reached at a 



The most recent updates may be 
|http: //www, ast . cam, ac .uk/^steurs/l 



found at 



mass of about 1000 A^q, and at a somewhat lower mass 
for stars with lower metallicity. 

If F = L/L^dd — > 1 in the interior of a 1-D stellar 
model, the radiativ e pressure can be balanced by a density 



inversion (see, e.g., iLanger 1997). In our calculations, the 



highest luminosity is reached for homogeneous models for 
M « 1001 M0. In this case, F 0.98 in the outermost 
meshpoint in the stellar model (see Fig. [1]). A further in- 
crease of the stellar mass is then not possible, since the 
convection in the surface layers is insufficient to transport 
the stellar luminosity, and density inversions cannot build 
up. 

If stars rotate and the critical velocity of rotation be- 
haves as Ucrit oc (1 — F)°-^, critical rotation is achieved 



before F reaches unity for any rotational velocity (jLanger 
1997). For rotating stars, the actual upper mass limit for 
solar metallicity may therefore be lower than 1000 M©. 

The development of a core-halo structure leads to a 
bending of the sequence of homogeneous stars towards the 
red in the HRD. In our models, this bending o ccurs above 
133 Mq in a g reement with e.g ., Figer et al 



a mass 



' ^.^^j^^^^ in a g reement with e.g . - 
(Il998h : llshn et all (Il999l ): IStothers fc ChinI (|l999^ . 



An approximate mass-luminosity relation for (100 - 
1000) Mq homogeneous stars obtained in the present 
study is (in solar units): 



^03.48^^1.34^ 



(1) 



For the (25 - 115) Mq initial mass range, the mass-radius 
relation may be approximated as 



R ^ 10-°-9m1°2, 
and for larger masses as 



R ~ 10 



-0.77^^0.96 



(2) 



(3) 



We note that although we can theoretically con- 
struct this sequence of completely homogeneous mas- 
sive star models, in reality, stars may form by accretion, 
and stars could start core hydrogen burning before ac- 
cretion halts, and evolve along the main-sequence un- 
til the accretion reservo ir becomes exhausted (see, e.g., 
Behrend fc Maedeill200ll ). 



2.2. The mass-loss rate 

Massive stars are subject to considerable mass loss, driven 
mainly by radiation pressure through spectral lines, which 
can be enhanced by currently poorly understood vibra- 
tional modes. For stars more massive than ^ 100 Mq 
theore tical models are not yet well developed. iKudritzki 
(|2002r ) calculated mass-loss rates for masses up to 3OOM0, 
but the mass-loss rates of Kudritzki do not include the 
important effect of multi-line scattering, which may be vi- 
tal fo r models in close proximity to the Eddington limit. 
Vind ( ,2006 ) presented mass-loss predictions for very mas- 
sive stars close to the Eddington limit, and found a steep 
behaviour of the mass-loss rate as a function of F (but for 
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a constant ratio of the terminal velocity over the escape 
velocity). 

We did not use the mass-loss rates of Vink et al. I (|200d) 
for stars above 100 Mq, as they were derived for F < 0.5. 
As the Vink et al. recipe includes bi-stability jumps, which 
are a function of TcS and wind density (and hence F), 
we did not extrapolate the Vink et al. recipe outside its 
validity range. 

A break of evolutionary sequences occurred in cal- 
culations for massive stars by other authors with mass- 
loss prescriptions based on radiation-driven w ind mod- 
els (e.g., for Pistol star bv iFiger et all [l998h . The ex- 
trapol ation of fits to the empirical data on mass-loss 
rates ( de Jager et al. 19881) for stars more massive than 



~ 120 M0 results in unreasonably high mass-loss rates, 
causing the stellar evolution models to lose their ability 
to converge. 

Due to the sparsity of theoretical expressions for the 
amount of mass loss from supermassive stars, we adopt 
the following ad-hoc "momentum" equation for our evo- 
lutionary models: 



M = 



1 



,c(l-F)("-0-5)^ 



(4) 



where L is the stellar luminosity, Vqo is terminal velocity 
of stellar wind, c is the speed of light, and F = L/LEdd 
in the outermost meshp oint of the stellar model. After 
Kudritzki fc Pulsl (|2000h we define 



C(Teff)[(2GM/i?)(l-F)] 



0.5 



(5) 



We do not specify a particular physical mechanism 
for our mass-loss prescriptions, but the luminosity depen- 
dence in Eq. Q accounts for the assumption that radia- 
tion pressure almost certainly plays a role. Furthermore, 
we account for the proximity to the Eddington limit by 
introducing the L/ Lp.dd dependence. We have tested our 
Eq. (Ill) against the lVink et al. ( 2000( ) mass-loss recipe for 
a 60 Mq star, and although there are differences due to 
the presence of bi-stability jumps in the Vink et al. recipe, 
which we do not account for in our mass-loss prescription, 
the overall mass lost (35 M©) is very similar, bringing a 
60 Mq star do wn to 25 Mp 

As noted by Vink et air ( 200lh . while the metal lines 



are responsible for driving the wind, hydrogen and helium 
lines are sparsely populated in the spectral range where 
early-type stars emit most of their energy. This provides 
some justification for using the same mass-loss recipe for 
hydrogen-rich stars, as well as for stars with hydrogen- 
exhausted atmospheres. 

We base our choice of a in Eq. (4) on two criteria: 
1) stars should spend most of their lifetime between the 
ZAMS and the HD-limit (as defined in Humphreys & 
Davidson 1994), and 2) they should spend only a limited 
amount of time, e.g. less than ^ 1 — 2% of their lifetime, 
at effective temperatures cooler than the Humphreys- 
Davidson limit. The nature of the Humphreys-Davidson 
limit is still elusive, as it is still unclear whether the lack 



of stars above this empirical limit is due to physical pro- 
cesses (likely to be related to the Eddington hmit), or it 
has a statistical nature instead. As some stars are observed 
on the cool side of the HD-limit (see Fig. [TT|l. we consider 
the possibility that the evolutionary tracks extend to the 
yellow and red regions of HR-diagram. 

We present our simulations for 60 and 120 Mq stars in 
Fig. [2] and Tabled 

For a = 0.25, the Mq = 60 Mq star loses 23.5 Mq 
during the core-hydrogen burning phase that lasts for 3.8 
Myr. At the TAMS, the surface hydrogen and helium 
abundances are virtually equal X RiO.49 and Y «0.49, 
and C/N and 0/N ratios are very close to their nuclear 
equilibrium values (see Fig. 1131 below). The star may be 
classified as a strongly helium-enriched 0-star. A further 
12 Mq are lost in the part of the track where the stellar 
effective temperature is below ~ 10 OOOK, with a lifetime 
in this later stage of 0.048 Myr. When the star crosses 
the HD-limit moving blueward, the hydrogen abundance 
at its surface is 0.096, and the object has become a WR 
star. 

This track agrees well wi th the track for wh ich we used 
the mass- loss prescription of lVink et al.l ( 2000l ). where the 
mass-loss behaviour is not completely continuous due to 
the presence of bi-stability jumps. The total amount of 
mass lost during the main-sequence phase is 19.4 Mq, 
whilst 13 Mq is lost when star is cooler than ^ 10 OOOK. 
The lifetime of the la ter stage is 0.088 Myr. Th is compares 
well to the results of lLimongi fc Chieffil (120071) . who used 
the mass-loss rates from Vink et al.l (|200o[ ) and found a 
total amount of mass loss of AM = 22 Mq during the 
main-sequence stage. 

The Mq = 120 Mq star, computed with a — 0.25, loses 
~ 57Mq over its 2.76 Myr long main-sequence lifetime. 
While on the cool side of the HD-hmit (at T^s ;S 30 OOOK), 
the star loses an additional ~ 14 Mq in 0.039 Myr. In 
the latter case, the star becomes a WR object during the 
main-sequence (see more below). 

The last models of the sequences for stars with ini- 
tial masses 60 and 120 Mq are in the core He exhaus- 
tion phase (Table O, with masses of ~ 40 Mq and 
~ 21Mq, respectively. These values are in reasonable 
agreement w ith the masses of 30 M q and ~ 20 Mq, 
obtained by iLimongi fc Chieffi _( 2007 ) who used mass- 
loss rates from Vink et al. I (l200fll )~ for blue supergiants, 
from de"jager et al.l (Il988l) for red supergiants, and from 



Nugis fc LamersI (|2000[ ) for WR stars. We a lso no te that 



the tracks computed by iLimongi fc Chieffil (|2007l ) for 60 
and 120 Mq stars penetrate into the cool region of the HR- 
diagram to log Toff ~ 3.65 and 3.70, respectively, again in 
good agreement with our computations. 

Finally, our tracks for a — 0.25 are in reasonable agree- 
ment with the track s computed using the spline fits of 
de Jager et al. ( 19881 ) for empirical M but, as we men- 



tioned above, the latter fits diverge at high stellar masses. 

For a = 0.1, the Mq = 120 A/q star crosses the 
Humphreys-Davidson limit during core-hydrogen burn- 
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ing, and spends ~ 0.5 Myr in the low-temperature region, 
which is too long. 

Based on this analysis, we conclude that a ~ 0.25 is 
suitable for evolving supermassive stars and we adopt this 
value throughout the paper. 

2.3. Evolution of supermassive stars 

A summary of the parameters of evolutionary tracks com- 
puted in the present paper is given in Table Figure [3] 
shows evolutionary tracks for stars with initial masses 60, 
120, 200, 500, and 1001 Mq, computed for different values 
of the parameter a. 

For a > 0.5, the mass-loss timescale is shorter than the 
nuclear evolution timescale. The stars evolve downward 
along the locus of homogeneous stars remaining practi- 
cally unevolved (Figs. I3l4p . For all sequences, the mass- 
loss rate is so high that nuclearly-processed layers are al- 
most immediately exposed at the surface. For stars with 
Mo=500 and 1001 M©, when central hydrogen becomes 
exhausted, the surface abundance of hydrogen is just sev- 
eral per cent. Thus, SMS turn into WR stars during the 
core hydrogen burning stage. 

For models with high a, the bending of the evolu- 
tionary tracks to higher Teff (bottom loop) starts when 
Xc — 0.05 in the convective core, i.e. this corresponds 
to what, in the common notation, would be "point B" of 
the evolutionary track. The blue points of the loops of 
the tracks for a > 0.5 correspond to the core hydrogen 
exhaustion. The direction of evolutionary tracks changes 
from redward to blueward, when the core He-burning be- 
comes the dominant source of energy: Lne/LH ~ 3. For 
all models, the stars evolve into virtually identical helium 
stars, as the stellar evolution converges. The evolutionary 
tracks start to deflect from the ZAMS to the right only 
for a < 0.5. These high-a tracks may be relevant for the 
evolution of putative SMS (M > 120 Mq), but they are 
less relevant for known VMS (M < 120 Mq). 

The behaviour of the evolutionary sequences w ith high 
mass- loss rates is consistent with the results of iMaeder 
(1983), who showed that stars with an initial mass in 
excess of 100 M©, and with a time averaged M ^ 2 x 
10~^ Mq yr~^ evolve towards the left of the HR diagram. 
In the interior, a large homogeneous core is surrounded by 
an envelope with very little difference in composition be 



tween the core and the stellar surface. IStothers fc Chin 
( 1999[ ) used OPAL opacities and found a similar conver- 
gence of the evolutionary tracks. 

To be consistent with restrictions upon the mass-loss 
rates found by evolutionary computations for Mq =60 
and 120 Mq stars (§ 12. 2p . we now discuss the tracks for 
a — 0.25. The behaviour of all evolutionary sequences 
is rather similar. For high a {a—l\ Fig. d]) the tracks 
converge quicker than for our favourite q;=0.25 models 
(Fig. [5]). Figure [5] represents the relations between evo- 
lutionary time, the mass of the star, the central hydro- 
gen abundance and the surface helium abundance for 



a = 0.25. For the Afo=200 and 1201/© models, the hy- 
drogen surface abundances for Xc = are 0.12 and 0.25, 
respectively. If we consider the surface helium abundance, 
Y w 0.4, as a conventional threshold for assigning the 
WR spectral type, this limit is reached within < IMyr for 
stars with M > 200 M©, whilst it takes about 2 Myr for 
the 120 M© star. 

Figures [6] and [7] show that after the end of the cen- 
tral H-burning all stars start to expand on the thermal 
time scale of their envelopes. Expansion is accompanied 
by a strong reduction of L/L^dd- although the stellar L 
remains practically the same, the LEdd of the outermost 
meshpoint increases rapidly due to a reduction in opacity. 
The reduction of F = L/L^dd would increase the mass- 
loss rate if it were to act alone. However, our adopted 
mass-loss prescription depends weakly on F. The decisive 
factor in Eq. ([4]) is i^oo , which decreases strongly with the 
expansion of the star, causing the growth of M. A further 
mass and gravity reduction at almost constant luminosity 
results in a kind of runaway mass loss, i.e., the star be- 
comes unstable. The expansion terminates when the core 
He-burning becomes the dominant source of stellar lumi- 
nosity. The star then returns rapidly to the blue region 
of the HR-diagram. In our calculations mass loss occurs 
continuously. In real stars, one may expect that instabil- 
ity associated with high L^dd would manifest itself as a 
sequence of sporadic mass-loss episodes which terminate 
when He-burning causes an overall contraction of the star. 
We may speculate that the total amount of mass lost in 
the "spike" of the mass-loss rate during redward excur- 
sions is similar to that which would be lost by a "real" 
star in a series of outbursts. The total stellar luminosity 
becomes very close to the Eddington luminosity when the 
star returns to the high Tofi region of the HR-diagram. 

An interesting feature of SMS evolution is the increase 
of the core hydrogen-burning time of the stars with in- 
creasing a in the mass-loss prescription. This is under- 
standable, as heavy mass loss is accompanied by a re- 
duction of the central temperature and density (like in 
close binary components that experience dynamical- or 
thermal-timescale mass loss upon Roche lobe overflow). 

Figure [8] shows that irrespective of initial stellar mass, 
for a — 0.25, the maximum mass-loss rates do not ex- 
ceed the maximum ma ss-loss rate for line-driven winds. 
Smith fc Owocki ( 2006h recently estimated this maximum 
to be: 



M « 1.4 X 10-*{L/{10^ Lq)) Mrs, yr~ 



(6) 



A similar order of magnitude (10"'^ Mq yr"-'^) estimate 
for the upper limit of M due to radiatioii pressu re on spec- 
tral lines was obtained by lAerts et all (|2004h . Figure [9] 
shows the run of the "wind performance number" 77 — 
{Mvoc)/{L/c), along evolutionary tracks for model se- 
quences with a = 0.25. 

One of the main results in the section is the finding of 
the convergence of our models with a wide range of ini- 
tial masses (see Fig. [5]). This convergence of stellar masses 
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which initially differ by a factor ^10 may crudely be ex- 
plained as follows. During most of the stellar lifetime the 
value of L/L^dd is not far from 1 (see Fig. [7]), whilst the 
stellar radius changes significantly for only a very short 
time. Therefore, the crucial factor that sets the rate of 
mass loss is the stellar luminosity, which may, as a first 
approximation, be taken as the luminosity of a homoge- 
neous star. For the case of a = 0.25, the stellar lifetime is 
t 10^-3 (M/Mq) -0-2 yr (for 120 Mq < M„ < 1000 Mq). 
When we now apply the mass-luminosity relation, the 
amount of mass lost, AM, is found to be comparable to 
the initial mass of the star, Mq. 

3. Models of VMS versus observed stars 

Figure [5] shows the mass-loss dependence on time for 
a = 0.25 sequences and compares them with the mass- 
loss ranges found for young massive stars in the Arches 
and Quintuplet clusters, and the cluster of Hel-emission 
stars in the Galactic center, as well as for a sample of 
WNLh stars in young massive clusters in the Milky Way 
and the Magellanic Clouds (A. de Koter, unpublished). In 
Figure [TOl we show in detail the relations between stel- 
lar mass and the mass-loss rate for a = 0.25 evolutionary 
sequences for objects with initial masses of 60, 120, and 
200 Mq , and we compare them with mass- loss prescrip- 
tions suggested for massive stars, as well as with obser- 
vations of known objects. We remind the reader that we 
consider stars with surface helium abundance Yc > 0.4 
and Tes > 10000 K to be WR stars. 

A fit for observed rn ass-loss rates of 0-stars is given 
bv Ide Jager et all (|l988l) . For M > 120 Mq, comparing 
models with any observational data fits is not very sen- 
sible, because it would involve too rough an extrapola- 
tion. However, our mass-loss rates for the Mq = 120 Mq 
star, when it is still H-rich, do not contradict the fits of 



de Jager et al. (|l988l) . nor the observations of the most 



massive H-rich WR-stars. 

For WR - starS; fits to observational data were given by 
e.g. iLangen l|l98£ 



3 



M = -2 X lO^^M^-^ 



Nelemans fc van den Heuvell (j200ir ): 

M = -1.4 X lO-^M^-*'^; 
de Ponder fc VanbeverenI (|2003lF I: 



log(-M) =logL- 10 



The numerical coefficient in Langer's formula is taken to 
be three times lower than the originall y suggested lowest co- 
efBcient, to account for wind clumping (|Hamann fc Koesterkd 
ll998l : lMarchenko et al.|[2007l ). 

* For the most massive W R stars, this fit gives ra tes close 
to the ones suggested by the iNugis fc Lamerg (|2000l ) fits. In 
Fig. 1101 the original formula is extrapolated to higher stellar 



(solar units and rates per yr are used). In Fig. \TU\ we 
also show observational estimates for mass-loss r ates o f 
WN stars with M > 20Mq from 'Hamann et al.| (|2006l ). 
Estimates of M for WR stars |Hamann ct_al. 2006 ) shown 



in Fig. [To] agree well with the estimates of iCappa et al 
()2004l) . based on radio-observations that showed that 
MwR hardly exceeds 10~^ Mq yr~^. Figure [10] clearly 
shows that extrapolation of the formulae given above to 
very massive WR stars is not justified. Late-type WN stars 
like WR20a (see Appendix for more details) may still be 
in their core-hydrogen burning stage. 

In Figs. [Tl] and [121 we present evolutionary tracks 
for a=0.25 and 0.5 sequences in the HR-diagram and we 
compare them to the positions of some of the most lumi- 
nous stars in the Milky Way and the Magellanic Clouds 
(Table [2]). For notes on the individual stars shown in 
Figs. [TT] and [121 we refer the reader to the Appendix. 
These figures indicate quite clearly that a = 0.5 would 
be too high a value, and that a should be ~ 0.25 to be 
consistent with the observations. 

The data presented in Figs. [6] [TO] [ill and [12] suggest 
that the initial masses of the stars observed in the Arches 
cluster. Quintuplet and R136 may exceed ~ 100 Mq. 
Shortly after leaving the ZAMS, these stars acquire sur- 
face abundances that result in their spectral classification 
as transitional types between 0-stars and WR stars. The 
positions of luminous stars with estimated M and M in 
Figs. [To] and [TTl suggest that they have initial masses up 
to 100 Mq, and are in the core-hydrogen burning stage. 
The positions of these objects in the (M — M)-plots and 
the HR-diagram are consistent. This point may be further 
explored when more detailed grids of evolutionary tracks 
become available. We note that masses and the estimates 



for M are both uncertain. For instance, iRepolust et al 



(j2004[ ) estimate the range of uncertainty of their spectro- 
scopic mass determinations for Galactic 0-stars (listed in 
Table [2]) as l3o%- For the mass-loss rates, the errors are 
up to 0.2dex. 

The mass-loss rates that we obtained for pure helium 
stars have the same range as M for hydrogen-rich stars. 
This is consistent with the similarity of M values found for 
0-stars and Hel-emission stars, of which s ome are helium 
riclo but it does not seem to agree with the lHamann et al. 
(|2006l ) data for hydrogen-poor WR stars. 

We should note that we obtain higher masses for WR 
stars and pre-supernovae masses than in some previous 
computations for non-rotating stars. Fo r instance, for the 
same 120 Mq star, ISchaller et all (|l992l ) obtain ~ 7.6 Mq 
at the end of their calculations, whilst our computations 
yield a final mass of 40.4 Mq (for a = 0.25). The difference 
stems from the difference in mass-loss rates for He-rich 
stars. For the latter, Schaller et al. use the rates from 
Langei uncorrected for clumping, that are much 

higher than those given by our empirical algorithm for M. 



^ The estimates of M for Hel stars strongly depend on an 
uncertain He/H-ratio. 



Yungelson et al.; Super-massive stars 



7 



In Fig. [131 we show the evolution of the chemical abun- 
dances at the surface of a 60 Mq VMS and a 500 Mq 
SMS. We note several differences in the behaviour of the 
surface abundances. On the surface of the Mq =&QMq 
star, helium starts to dominate over hydrogen only close 
to the TAMS, whilst for the Mq =500 Mq SMS, this al- 
ready happens much earlier. This difference is due to the 
relative sizes of the convective cores. Furthermore, the sur- 
face of the VMS becomes enriched in nitrogen (relative to 
the solar abundance) close to the TAMS, and remains so 
until the end of the helium exhaustion of the core. For this 
case, we expect the pre-supernova object to be a WN star. 

The surface layers of the SMS become enriched in ni- 
trogen almost immediately after the departure from the 
ZAMS. Later, at the TAM^, the hydrogen and nitrogen 
abundances drop almost simultaneously, and the surface 
becomes dominated by He and C. In the final stages of 
evolution, oxygen becomes the dominant surface element, 
whilst the abundance of helium becomes even lower than 
that of carbon. In this case, the expected pre-supernova 
object is an oxygen-carbon star, with traces of helium at 
the surface. 

Figure [13] also shows the differences in the surface 
CNO-cycle element ratios of VMS and SMS. For the lat- 
ter, the layers in which the elements are in nuclear equilib- 
rium are exposed virtually immediately at the beginning 
of their evolution (after the loss of only about 100 Mq in 
only ~0.7 Myr). This situation persists until C and O en- 
riched layers become exposed due to He-burning. As the 
core H-burning evolutionary stage is by far the longest, 
we conclude that a majority of SMS should have nuclear- 
equilibrium ratios of CNO-elements. 

For the case of the VMS, the transition in chemical 
abundance from solar to nuclear-equilibrium occurs some- 
what more smoothly. It starts after about 2 Myr of evo- 
lution, whilst it takes about 2 Myrs out of a total 6 Myrs 
MS-lifetime. During this transitional phase, the star loses 
about 10 Mq. If our mass- loss algorithm is applicable to 
the lower mass VMS, we should expect this transition to 
occur smoothly. VMS close to the TAMS are expected 
to be associated with luminous blue variables (LBVs), for 
which a w ide span of CNO-el ements ratios has been found. 
Although ls"mith et al.l (| 19981 ) found LBV ejecta to be N- 
enriched, these have generally not yet reached CNO equi- 
librium values. However, as there is controversy about the 
origin of LBV nebulae, it would arguably be better, or 
at least more direct, to consider photospheric abundances 
instead. These, however, may depend on complexities of 
the atomic physics and wind analyses, and the results are 
conflicting. For instanc e, evidence for advanced nuclear 
processing was found bv lLennon et al. ( 1993h for the LBV 
R71 in the LMC, but this is not a well-established fact for 
LBVs as a population. Although most LBVs show CNO- 



enriched atmospheres, it seems unlikely that all of them 
have already reached their CNO equilibrium values. 

A range of CNO abundances may naturally be ex- 
plained from our VMS computations, but given that most 
LBV masses are not significantly over 150Mo, it seems 
unlikely that all LBVs are the progeny of SMS, as our 
models would predict SMS to reach CNO equilibrium al- 
most instantaneously. 

4. The fate of supermassive stars 

According to our evolutionary calculations, SMS trans- 
form into oxygen-neon stars with traces of carbon. 
Depending on the progenitor mass, the final masses of the 
stars are between ~20 Mq and ~140AfQ. One may ex- 
pect that further reduction of the total stellar mass would 
be negligible, because of the short lifetime of stars in their 
post core-helium burning evolutionary stages. The helium- 
carbon-oxygen outer layer comprises only a few per cent 
of the total mass of the star. Thus, one may expect that if 
these objects were to explode, SMS would probably pro- 
duce type Ic events. 

It ha s been known for decade s (see, e.g., pioneering pa- 
pers of iRakavv fc Shavivlll967t lFralevlll968l : lOber et al 



With the caveat that we use the definition of the TAMS in 
a broad sense indicating the end of the core hydrogen-burning 
stage. 



19831 and many other studies) that if the mass of oxygen 
stars exceeds several tens of Mq, during central oxygen- 
burning, they enter the electron-positron pair instability 
regime, and contract quasi-dynamically. When the central 
temperature increases to (3 — 6) x 10^ K, central oxy- 
gen burning becomes explosive, which is much faster than 
neutrino energy losses. The released nuclear energy may 
be sufficient for the internal energy to exceed the gravi- 
tational binding energy. The star will then disrupt com- 
pletely, without leaving a compact remnant, giving rise to 
a so-called "pair-instability (or pair-creation) supernova" 
(PISN). If the released energy is not large enough to dis- 
rupt the star, the star collapses to a black hole. This insta- 
bility arises irrespective of the metallicity of the progenitor 
star. However, for solar composition stars the possibility 
of PISN is usually not considered, since in the evolution- 
ary models for even the most massive stars formed in the 
"standard" way (i.e. without possible runaway collisions), 
with masses up to "observational" upper limit of 150Mq, 
sufficiently massive oxygen cores are not formed. We note 
that this may stem from an overestimate in the mass-loss 
rates during the WR-stage (if mass-loss prescriptions like 
those shown in Fig. [TO] are extrapolated). 

However, massive oxygen cores are more often consid- 
ered in relation to the evolution of Population III stars. 
Massive oxygen cores may experience only moderate mass 
loss, owing to the absence of metals (but see Vink & de 
Koter 2005 and Vink 2006 for a discussion on the pos- 
sibility of higher mass loss due to the radiative driving 
through intermediate mass CNO elements, and the likely 
proximity to the Omega/Eddington limit). One may use 
the results obtained for Pop. HI stars for solar composi- 
tion stars, because of a rather modest difference in inter- 
nal structure of massive stars of different initial metallicity 
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during the late stages of evolut ion (e. g., ISchaller et al 
19921) . lUmeda fc Nomotol (|2002l ) find, in computations for 
Pop. Ill stars, that PISNs are experienced by stars that 
form He-cores with mass (70 - 129) M(7). A similar range 
of (63 - 133) Mq was found by iHeger fc WooslevI (|2002h 
from computations of non-rotating He-stars. Apart from 
differences in the input physics, some discrepancy in mass 
obtained by different authors stems from the well-known 
fact that the evolution of initially naked helium stars is 
slightly different from the evolution of similar mass helium 
cores of stars formed through hydrogen burning. Stars 
that have helium cores more massive than the ranges of 
PISN-progenitors quoted above form black holes without 
ejecting matter. Stars with He-core masses below these 
limits are expected to collapse direc tly to a black hole or 
form a black hole through fall-back ( Frver| [l999). 

In Fig. [HI we plot the initial-final mass relation 
for SMS, using the re s ults of lUmeda fc Nomotol (|2002h 
and iHeger fc WooslevI (j2002l ). and we mark approxi- 
mate boundaries between regions with different outcomes. 
These boundaries are rather crude, but they suggest the 
following: we hardly expect the formation of black holes 
with masses larger than ~ 150 Mq. As the difference be- 
tween the upper mass boundary of PISN-producing SMS 
(~ 800 Mq) and the upper SMS initial mass hmit 
1000 M0) is quite small, black holes with masses larger 
than PISN-producing objects may hardly form, and it is 
quite likely that the most massive black holes produced 
by SMS have only M < TOMq. 

Figure [15] shows the lifetimes of supermassive stars. 
It is worth noting that lifetimes of SMS may be shorter 
than the 3 Myr usually assigned to the most massive stars, 
and this may be relevant for studies of the upper li mit of 
nasce nt stars that use age-related arguments (e.g. 
200i). 



Figei 



5. Summary and discussion 

In this paper we discuss the evolution of stars with initial 
masses in the range 60 Mq to 1001 M©. Our study was 
motivated by the results of the direct iV-body simulations 
of dense star clusters by Portegies Zwart et al. (1999), in 
which a star grows by repeated collisions to well beyond 
100 M0. In later studies, Portegies Zwart et al. (2004) 
proposed that the collision runaway in star clusters could 
explain the presence of a black hole of ^ 600 Mq in the 
star ULX M82 X-1 in cluster MGGll, which supposedly 
could have formed from a ^ 1000 Mq star. According to 
these models, the mass of the VMS increases over the stel- 
lar lifetime, starting as a homogeneous massive ~ 100 Mq 
star that grows to Xl, 1000 Mq within the core-hydrogen 
burning stage of evolution. 

We have assumed that our stellar evolution calcula- 
tions are representative for stars that grow in mass via the 
collision runaway process. This is not necessarily correct, 
as the hydrogen reservoir of the collision runaway prod- 
uct will continuously be replenished by repeated collisions, 
whereas in our simulations we start with a high-mass ho- 



mogeneous model. Furthermore, rotation may play a rele- 
vant role in the evolution of these objects, which was also 
ignored in our study. Having noted this, we may nonethe- 
less expect these massive objects to be subject to heavy 
mass loss, and we may provide meaningful results in terms 
of the fate of the objects under consideration. 

We have confirmed previous results on the existence 
of an upper mass limit for chemically homogeneous stars, 
which is reached when the luminosity in the outermost 
layers of the stars approaches the Eddington luminosity, 
at which point gravity is unable to balance radiation pres- 
sure. For non-rotating solar composition stars, this limit is 
reached at about 1000 Mq. We evolved the stars of 60 Mq 
to 1001 Mq to the end of the core helium-burning stage 
and calibrated the stellar mass-loss prescription adopted 
by enforcing the condition that no star should spend more 
than a few percent of its life above the HD- limit. We found 
that mass-loss rates and HR-diagram positions of the ear- 
liest 0-type stars, and stars classified as transitional be- 
tween O- and WR-stars, may be consistent with our com- 
puted tracks. 

Based on our mass-loss prescription and stellar evo- 
lution calculations, we argue that the observed massive 
stars in the Galaxy and the Magellanic Clouds had birth 
masses of up to ~ 200 Mq. For stars in the lower part 
of this mass-range, mass loss during the MS-stage leads 
to the exposure of layers moderately enriched in He, and 
they may be observed as late-type WN stars with hydro- 
gen features in their spectra. In the upper part of this 
range, stars become hydrogen-deficient WR stars already 
on the main- sequence. 

Recentlv. iBelkus et all ( 2007 ) published a study of the 
evolution of stars with masses up to 1000 A/q. Our study 
differs from that of Belkus et al. in two important aspects: 
(i) they provided a simple evolutionary recipe based on 
similarity theory, assuming that stars are homogeneous 
throughout the entire course of their evolution due to vig- 
orous convection and stay in thermal equilibrium, whereas 
we present detailed stellar structure models; (ii) Belkus et 
al. used extrapolated mass-loss rates as predicted from 
radiation-driven wind theory, whilst we employed a mass- 
loss prescription that is consistent with the location of the 
most massive stars in the Hertzsprung- Russell diagram. In 
both studies, it is found that SMS are subject to dramatic 
mass loss that probably inhibits the formation of IMBH 
by the runaway stellar collision scenario. 

Star clusters that experience core collapse before the 
most massive stars have left the main-sequence can de- 
velop a supermassive star via collision runaway. The mass 
of such an object is accumulated in subsequent collisions 
on a time scale of less than 3 Myr. The mass which can be 
grown in this time interval can be estimated using Eq. (2) 
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of iPortegies Zwart et al] (|2006t fl: 



0.01m 1 



tr 



100 Myr 



-1/2 



(7) 



where mm is the mass of the object formed by run- 
away colhsion, m is the system mass, and tri is the 
system relaxation time. The average mass increase 
per c ollision is about ~ 20AfQ (.Portegies Zwart et al . 
I999I) ffor a more complet e disc ussion, however, see 
Portegies Zwart fc McMillanI (|2002l )). The supermassive 
star that accumulates ~ IOOOM0 has experienced some 
45 collisions between the moment of gravothermal collapse 
of the cluster and the moment that the supermassive star 
dies. The mean time between collisions for this model is 
^ 6.0 X 10^ years, resulting in a mass accretion rate of 
^ 3 X lO~^M0/yr. This rate of mass accretion is lower 
than the rate of mass loss that we get for the most mas- 
sive stars 3.8 x 10"^ Afo/yr), but since these numbers 
of mass-accretion and mass-loss rate are quite uncertain, 
it is not inconceivable that, in spite of the strong stellar 
winds, the objects may nonetheless experience a net gain 
in mass during their li fetimes. This conclusion was drawn 
bv lSuzuki et al.l ( 2007 ) who recently studied stellar evolu- 
tion with mass loss of collisionally merged massive stars. 
They concluded that stellar winds would not inhibit the 
formation of SMS. We note, however, that the Suzuki et 
al. (2007) calculations were limited to masses up to ^ 
100 Mq and their calculations did not consider the more 
important effect of mass loss for the final mass growth 
process, where the mass is supposed to increase beyond 
this value of 100 Mq by an order of magnitude. 

In this study, we have found that a super-massive star 
is likely to shed most of its mass well before it experiences 
a supernova explosion. In several cases the supernova pro- 
genitor still had ^ 100 A/q, and such objects could possi- 
bly collapse to black holes of intermediate mass but with 
Mbh « 1000 Af0. Our calculations suggest the possibility 
of the formation of objects that experience pair-instability 
supernovae, which would be interesting phenomena to 
observe. The recently discovered, ver y luminous super- 
nova 2006gy may have been a PIS N (jsmith et all 12007 : 
Ofek et al.ll2007l : iLanger et al.ll2007^ ■ although this is stiU 
under debate. 

Obviously there is the possible caveat that the quanti- 
tative difference between starting as a homogeneous high- 
mass star (as we considered in our paper), or growing 
to one over the main-sequence lifetime by repeated col- 
lisions may be significant, but at this stage we cannot 
provide further insights about the consequences of these 
potential differences. We therefore draw our conclusions 
on the calculations at hand, and we argue that the major- 
ity of supermassive stars probably end up as black holes 
of Af < 7OAf0, with the possible exception of some of 
them exploding as pair-instability supernovae. With the 



These results were obtained using standard models for run- 
away mass accumulation without accounting for the possibility 
of vigorous mass loss as considered in the present paper. 



approximate nature of our applied mass-loss algorithm in 
mind, an d also taking into acc ount the (also approximate) 
results of iBelkus et al.l (|2007l ). we infer that the accuracy 
on the limit for black hole masses is ~ ±30 Af0 . We there- 
fore conclude that most supermassive stars end their lives 
as 70 ± 30 Af0 black holes. 
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APPENDIX: Notes on individual stars 

We first discuss the enigmatic LBV 77 Car. We adopt a lu- 
minosity of 5 X 10^ Lq, which corresponds to the observed 
IR-flu x and assumed distance of 2.3 Kpc (Hillier et al 



l200lh . The positions of the two labels for the star in the 
HR-diagram reflect the ambiguity in the determination 
of the stellar radius, caused by the presence of an opti- 
cally thick wind that dominates the spectrum and pre- 
vents the determination of Ri,. The maximum and min- 
imum Toff correspond to a Rosseland optical depth of 
T — 155 and 0.67 respectively in t he atmosphere. 77 Car 



may be a 2020 ± 5 day binary (jPamineli et al.l 11997 



I2OOOI ). Assuming that the luminosity is Eddington, Hillier 
et al. (2001) estimate the minimum total mass of the 
system to be 150[L/(5 x 10^)] Mq), but some contribu- 
tion of the secondary star is probably not very signifi- 
cant, although the secondary is possibly a late-O/WR star 
(|Iping et al. I l2005h . Our computations show that the lu- 
minosity of a several hundred solar mass star may still 
be far from the Eddington limit. The estimated mass- 
loss rate is Af « 10"'^ M0 yr^-'^. Bearing in mind that 
A/q = (200 — 300) Mq stars in their core-hydrogen burning 
stage evolve at almost constant luminosity, L = 5 x 10^ Lq 
in our co mputations correspo nds to a ZAMS mass of about 
250^^0. iHillier et all (|200l[ ) derive for the rj Car atmo- 
sphere a H/He fraction of «5 by number, i.e. X«0.56, 
YwO.44. Thus, i, X, Y, M, and the "high" Toff solution 
suggest that 77 Car is a ^ 250 Af0 star somewhere in the 
middle of the core-hydrogen burning stage. The nature of 
its instability has yet to be found. 
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The position of the candidate LBV Pistol s tar is drawn 
accord ing to the "low-luminosity" solution of iFiger et al. 
(|l998l ). since iNaiarro et al. bund that its metal 



content equals Z ~ SZ©, corresponding to its location in 
the star-forming region of the Galactic center, and that 
this favours the lower luminosity solution for its stellar 
parameters. Note however that, if it turns out that the 
Pistol star is an LBV, its To// will be variable. Of all 
our tracks, the one for Afg = 200 Mq se ems to prov i de the 
best fit of the position of the Pistol star. lFiger et al. I (ll998l) 
find Y=0.3-0.4 at the surface of the Pistol star, as well as 
strong N enhancement. Going red ward, at log Toff = 4.15, 
our Mq — 200 evolutionary sequence has a mass 
M = 92 Mq and Y=0.86. Going blueward, it has a mass 
of M=86.5M0 and Y=0.91. We therefore cannot match 
the chemical composition of the Pistol star. Other tracks 
passing through the same position also have too high a he- 
lium fraction Y. For instance, the track for Mq = 120 AIq, 
which runs slightly lower than the error box of the Pistol 
star, has Y> 0.7. Note that the determination of the He- 
abundance is challenging. However, if the high metallicity 
of t he Pistol star is r eal, it may be an almost unevolved 
star: Ishii et al. ( 19991 . see their Fig. 6) find that the Pistol 
star lies very close to the ZAMS for stars with X=0.7, 
Z=0.1; its mass is then slightly higher than 200 MgH. 

There is another LBV star in the Quintuplet clus- 
ter, similar to the P istol star: FMM362 (jPiger et al.lll999l: 
iGeballe et aTl I2OOOI ). The latter authors report T^s = 
(10000 - 13000) K for this object, with L > 10^ Lq. 

For the stars in the Hel cl uster in the Galactic cen- 
tre, we show M estimates from lMartins et aL ( 2007 ). The 
positions of the Hel stars in the HR-diagram are consis- 
tent with their interpretation as being evolved blue super- 
giants, close to the WR stage. We note that the evolution- 
ary stage of the Hel star has to be very short. 

We note that the WN stars - termed AdK - may be in 
their core-hydrogen burning stage, as they are positioned 
along portions of the tracks where Y is only moderately 
enhanced. 

The cool hypergiant IRC +10 420 is the only object 
that is believed to be observed in the phase of rapid tran- 
siti qn from the red supergiant stage to the WR p hase (see 



e.g. iHumphrevs et al 



1997; Blocker et al. 19991 and ref 



ttpr 

m 



erences therein). Its T^g appears to have increased by 
1000 - 2000K wi thin the last 20 years. We note that 
Smith et al.l ( 20041 ) computed radiation-driven wind mod- 



els for objects with a low stellar mass and discussed the 
possibility that these cool hypergiants like IRC -flO 420 
may in fact be "LBVs in disguise", where a large mass- 
loss rate induced by the high Eddington factor produces 
a pseudo-photosphere. Although the object is located in 
the HR-diagram at the lowest luminosity level, based on 
our calculations we expect the existence of similar more 
luminous stars, which have to be extremely rare due to 



We must note however that the models o f llshii et all (|l999l ) 
for mass higher than ~ 200 Mq have more extended envelopes 
than our models. 



their short lifetimes, p Gas may be an example of star 
that is in an unstable post-RSG stage of evolution, when 
it experiences o utbursts ejecting mass at a rate of several 
0.01 Mq yr-i (|Lobel et al.ll2003l ). 

We show the position of one of the bri ghter 
LBVs, AG Car, after iHumphreys fc DavidsonI (|l994 ) and 



Lamers et al.l (|l996al ) . Two data labels illustrate the range 



of excursions over the HR-diagram experienced by the 
LBV. 

Especially re markable is WR20a, a Pprb ~ 3.68 day 



mass ive binary (jRauw et al.l |2004 l2005t iBonanos et al 



I2OO4I ) with both components of a WN6ha spectral type 
(see Fig. [TT|) . The compone nt mass estimates ar e 83 ± 



5.0 Mq and 82 ± 5.0 ( Bonanos et al.l 12004 ). This 



makes these components of WR20a the most massiv e 
stars weighed in binaries. According to lRauw et al. I (l2005l) . 
fundamental parameters of each of the components are: 
Toff = 43000 ± 2000K, Tboi/T© ^ (1-15 ± 0.15) x 10^ 
M — 8.5 X 10~^ Mq yr~^ (assuming a clumped wind with 
a volume filling factor / = 0.1). Nitrogen is enhanced in 
the atmospheres, whilst carbon is depleted. Spectral clas- 
sification implies only a weak helium enrichment of the 
atmosphere. The ori gin and evolution of this system de- 
serves urgent study. iRauw et"al ( 20051 ) propose that the 
position of the WR20a components in the HR-diagram 
"suggests that they are core hydrogen burning stars in 
a pre-LBV stage and their current atmospheric chemical 
composition probably results from rotational mixing that 
might be enhanced in a close binary compared to a sin- 
gle star of the same age" . However, we should note that 
the evolution with complete mixing of a high-rotational 
velocity star in a synchronised binary with mass loss may 
result in evolutionary tracks that evolve to the left of the 
ZAMS, in apparent conflict with the position of WR20a. 
As Figs. Ho] and [n] show, Toff, L, and M of WR20a are 
actually quite consistent with a star with an initial mass 
of about 100 Mq in a rather early stage of core-hydrogen 
burning, when the layers enriched in He are beginning to 
be exposed. 

In addition to the Galactic stars, we plot some 
Magellanic Clouds stars in Fig. [TlJ Most remarkable are 
the R136 cluster members of 03(If*) and 03If*/WN6A 
subtypes that are consid ered to be transitional b etween 
O and WN type stars (|Massev fc Hunted \l99^ . The 
stars are located along the ZAMS or slightly to the left 
of it. The R136 cluster belongs to the LMC and has a 
several times low er Z than Galactic stars. Figure 6 of 
Ishii et al.l (|l999f) shows that in the (100-200) M© range, 
ZAMS stars with Z=0.004 are 0.1 dex in logTes hotter 
than Z=0.Q2 stars, but they have the same luminosity. 
As our Fig. [TT] shows. (100-200) Mq stars may easily re- 
veal over-abundances of helium, whilst they evolve over 
AlogTcff ~ 0.1 in the HR-diagram. The surface helium 



^ The degeneracy with respect to the efltective tempera- 
ture is caused by the absence of individual determinations of 
Teff. The latter we re assigned according to spectral type after 
IVacca et all (|l996l ). 
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abundance increases to F = 0.4 in ^ 0.5 My r. Therefore, 
R136 cluster members may be very massive stars that 
are still in their core H-burning phase, with a strong 
stellar wind that causes the spectrum to mimic that of 
evolved WR stars, as was originall y conjectured fo r one o f 
these stars (R136-006=R136a3) bv lde Koter et all (|l997t ). 
Another possibility is that the surface layers enriched in 
helium are already exposed (this may happen in ~ 1 Myr, 
see Fig. [5]). This is also consistent with the star formation 
histo ry of R136 (see the discussion in Massev fc Hunter] 
19981). We also note t hat the effective-temperature scale by 
Martins et al. (20o3) assigns 031-type stars logTe// ~ 4.6 
(instead of.4.7 in Vacca et al. 1996) and smaller bolomet- 
ric corrections (implying a reduction of the luminosity by 
~ 0.25 dex). This would move these stars to the right 
of the ZAMS, in the proper direction. The position of 
R136 members in the HR-diagram is roughly consistent 
with the estimate of the upper limit of mas ses of stars in 
this complex of (140-160) Mq obtained bv iKoenI (|2006h 

by power-law distribution fitting to these stars. 

Finally, we also plot P Cygni ( Lamers et al. 1996b[ ) in 
Figs. [Tl] and [121 and show the positions of several other 
high-luminosity Galactic stars for which our calculations 
may be relevant (see Table [2] to identify them if they are 
not annotated in the Figures). 
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Table 1. Summary of calculated evolutionary tracks. Successive columns give the initial mass A/q, the value of the 
a parameter, the age at the terminal age main-sequence txAMS, the mass at the TAMS MtamSi the age of the last 
computed model tf, the mass of the last computed model Aft, the central abundances of He, C, and O in the last 
model Ycf, Xc,ci, -''^o.cf , and the abundances of H and He at the surface of the last computed model Xgf, Ygf. Note that 
helium burning results in the formation of oxygen-neon cores with trace s of carbon . The "V " in the second column 
stands for the 60 Mq track computed using the mass-loss prescription of IVink et al. I (|2000l) . 



iVIo 


en 


CTAMS , J^U yr 


Mtams 


tf,lU yr 


IVli 


-icf 


^C,cf 






V" ^ 
-i sf 




1 n 

i.U 




zy. DO 


KQA 

z.oy4 


op; 0/1 

Z0.Z4 


U.D4o 


n 0Q7 
U.Zo ( 


n HQQ 
u.uoy 


n n 

U.U 


n Q7P. 
u.y ( D 


1 nni 


u.y 


z.odU 


o4.oz 


71 Q 
Z. i io 


on op; 
zu.zo 


n ni Q 
U.Uio 


n 1 /1 1 

U.i4i 


n Qi 1 
U.oii 


n n 
U.U 


n oo/i 

U.zU4 


1 nni 


U.o 


OA A 


4i.4D 


^A^ 
Z.04i 


op; 1 1 
ZO. ii 


n nfio 
U.Uoz 


n OAO 
U.ZDZ 


n fiQn 
U.DoU 


n n 
U.U 


n oQi^ 

U.ZOD 


1 nni 


u.o 


i.yoo 


ftfi 81 
oD.oi 


01 s 
Z.Zio 


/I Q 1 7 

4y. i / 


n noQp; 
u.uuoo 


n (^d^ 

U.UDi 


n Qp;f^ 

U.oOD 


n n 

U.U 


n 1 7/1 

U.i ( 4 


1 nni 


U.o 


1 VQO 
i. iifZ 


01 O O 

ziy.y 


1 70Q 

i. / yo 


01 Q /I 
Zio.4 


n OQ 1 
u.yoi 


n nnnoQ 
U.UUUzo 


n nnnoo 
U.UUUzZ 


n nop; 
U.UzO 


n op;q 
u.yoo 


1 nni 


U.zO 


1 777 
I. i 1 1 


ooo n 
zzy.u 


ni o 
z.uiy 


1 QO fi 

ioy.D 


Q A 1 n~4 
y.4 • iu 


n n 
U.U 


n Qi o 

u.oiy 


n nQp; 
U.UoO 


n onn 
u.zuu 




1 n 

i.U 


9 p;i Q 

z.OiO 


OQ Kd 

zy.oD 


fiQQ 

z.oyy 


1 Q fil 
iO.Di 


n n 
u.u 


n 1 1 p; 

U.iiO 


n Q/i/1 

U.o44 


n n 

U.U 


n OQO 
u.uou 


oyo 


1 n 
i.U 


z.O / z 


oo on 
zy.zu 


p;oi 
z.oyi 


07 Q1 

z ; .yi 


n op;q 
u.yoo 


n noi 
U.Uzi 


n nno/i 
U.UUU4 


n nQQ 
U.Uoo 


n 0/1 Q 
u.y4o 


p;nn 
oUU 


1 n 
i.U 


z.oUO 


OQ Qd 

zo.yo 


AOQ 
Z.DZo 


zD.yo 


n a A Q 
u.y4o 


n no/1 
U.Uz4 


n nno/i p; 
U.UUU40 


n nQQ 
U.Uoo 


n 0/1 Q 
u.y4o 


pinn 

OUU 


n Q 
u.y 


Z.4oO 


'i'i Q1 

oo.yi 


aQ7 

Z.oO ( 


on nQ 
zu.uy 


n ni 
u.uiy 


n 1 /I p; 

U.i40 


n QO^i 

U.oUD 


n n 

U.U 


n 01 Q 
u.zio 


p;nn 
oUU 


U.o 


Z.O i Z 


Ad no 
4U.Uz 


AOA 


OQ 


n CiAfi 

U.U4D 


n 1 op; 
u.iyo 


n 7Q1 
U. ( oi 


n n/in 
U.U4U 


n OQO 
u.zoy 


p;nn 
oUU 


U.O 


1 n/1 

Z. iU4 


(D. iK) 


QOQ 


/I op; 
4D.yo 


n n 
U.U 


n n^^o 
U.Udz 


n Q p;o 
U.oOz 


n n 
U.U 


n 1 op; 
U.iyo 


pinn 

OUU 


u.o 


1 Qfil 

i.yoi 


IDO. ( 


1 QSO 

i.yoz 


iDO.O 


n QQ1 
u.yoi 


U.UUU40 


n nnnoQ 

U.UUUZO 


n np;7 
u.uo / 


n ooQ 
u.yzo 


500 


0.25 


1.956 


183.4 


2.206 


115.1 


0.41 ■ 10~^ 


0.0 


0.827 


0.039 


0.215 


200 


1.0 


2.951 


26.76 


2.997 


25.74 


0.960 


0.014 


0.00032 


0.051 


0.930 


200 


0.9 


2.837 


30.75 


3.208 


18.98 


0.010 


0.128 


0.828 


0.0 


0.228 


200 


0.8 


2.740 


35.51 


3.083 


22.31 


0.027 


0.157 


0.786 


0.0 


0.238 


200 


0.5 


2.511 


56.71 


2.525 


54.61 


0.981 


0.0004 


0.00024 


0.047 


0.934 


200 


0.3 


2.379 


81.92 


2.397 


76.12 


0.980 


0.0004 


0.00024 


0.073 


0.907 


200 


0.25 


2.350 


92.28 


2.629 


61.88 


0.0 


0.00556 


0.846 


0.0 


0.304 


120 


1.0 


3.345 


24.15 


3.360 


23.67 


0.974 


0.0012 


0.00027 


0.074 


0.398 


120 


0.5 


2.906 


44.36 


2.915 


44.19 


0.975 


0.00070 


0.00025 


0.094 


0.887 


120 


0.4 


2.836 


50.74 


2.844 


50.52 


0.980 


0.00037 


0.00025 


0.133 


0.848 


120 


0.35 


2.813 


54.41 


3.108 


37.86 


0.981 


0.00036 


0.00025 


0.163 


0.818 


120 


0.3 


2.771 


58.79 


3.080 


39.34 


0.847 


0.073 


0.847 


0.0 


0.320 


120 


0.25 


2.786 


63.02 


3.077 


40.38 


0.043 


0.146 


0.781 


0.0 


0.427 


120 


0.10 


2.900 


47.64 


3.007 


44.48 


0.010 


0.077 


0.851 


0.077 


0.548 


60 


0.5 


3.394 


29.56 


4.357 


18.34 


0.0 


0.123 


0.834 


0.0 


0.975 


60 


0.25 


3.807 


36.31 


4.213 


20.94 


0.0 


0.113 


0.834 


0.0 


0.978 


60 


V 


3.685 


40.60 


3.787 


26.20 


0.719 


0.228 


0.023 


0.0 


0.876 
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Table 2. Some of the most luminous stars in the Galaxy and Magellanic Clouds. 



15 



No. 


Star 


Sp. type 


logL/L© 


logTeff 


M, IQ-'^ Mq yr"^ 


M/Mq 


Ref 


Comment 


1 
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T MP 
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4. DO 
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4K 


Q 





1 A 
r^ 
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4.D40 
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1 p; 
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ijiiiui. vvo-iy 


Uz il 


fi Q/1 


A fi/1 Q 

4.D40 


on 
ZU. 


1 QQ 

lyo. 


Q 





16 


K13D-036 


Uz It 


5.7 


4.633 


14. 


31. 


3 




17 


K13D-020 


02 If* 


5.9 


4.628 


23. 


An 

40. 


3 




18 


Sk -67 22 


02 If* 


5.69 


4.623 


15. 


23. 


3 




19 


LH90:Br58 


03If/WN6 


5.9 


4.613 


40. 


40. 


3 




20 


R136-014 


03.5 If* 


5.9 


4.58 


23. 


53. 


3 




21 


Sk -65 47 


04 If 


5.97 


4.60 


12. 


61. 


3 




22 


AV 75 


05.5 1(f) 


6.16 


4.60 


3.5 


96. 


3 


SMC 


23 


AV 26 


06 1(f) 


6.14 


4.58 


2.5 


91. 


3 




24 


HD93250 


03 V ((f)) 


6.01 


4.662 


3.45 


83.3 


4 




25 


HD66811 


04 1(f) 


5.90 


4.613 


8.8 


53.9 


4 




26 


HD 14947 


05 If+ 


5.90 


4.574 


8.52 


30.7 


4 




27 


HD15558 


05 Ill(f) 


5.93 


4.613 


5.58 


78.7 


4 




28 


HD210839 


06 I(n) fp 


5.83 


4.556 


6.85 


62.2 


4 




oo 


xliJoUD14 




0.00 


A /1fi9 
4.4DZ 


D.U4 


Q7 fi 

0/ .D 


4 




■^n 

ou 
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iVi vv 


31 


Rn+fiO 261 3 


M3 I 




3.544 






5 




32 


Mk51 


03I*f/WN7-A 


6.23 


4.648 






5 


LMC 


33 


R139 


07Iafp 


6^40 


4^558 






6 




34 


Mk26 


04 Ill(f) 


6.02 


4.647 






6 




35 


Mk24 


03V 


6.00 


4.686 






6 




36 


Mkl4 


03-6V 


5.94 


4.646 






6 




37 


IRC-I- 10420 


OH/IR 


5.7 


3.845 


300-600 




7 


MW 


38 


CygOB2 No.9 


05f 


6.6 


4.650 


12.7 


160. 


8 




39 


CygOB2 N0.8A 


05.5 1(f) 


6.18 


4.585 


13.5 


118. 


8 




40 


CygOB2 No. 12 


B8 la 


6.20 


4.049 


38.5 


71. 


8 




41 


HD 33579 


A3 Ia+ 


5.72 


3.902 


2. 


20.-30. 


9 


LMC 


42 


HD80077 


B2 Ia+ 


6.4 


4.23 






10 


MW, (LBV?) 


43 


HD119796 


G8 Ia+ 


5.7 


3.67 






10 




44 


HD152234 


B05 Ia(N wk) 


5.87 


4.41 


2.7 




11 




45 


HD15236 


B1.5 la 


6.1 


4.26 


6.0 




11 




46 


HD224914 (p Gas) 


F8 Ia+ 


5.70 


3.86 


54000. 




11, 12 


. . . , M in outburst 



Refe rences: 1 - iKo jjnigsUergerj |2004l). 2 - iMassev^eFaLl (Il989fl. 3 - iMassev et all (120051 1. 4 - iRepolust et al.l ( 2004). 
5 - iLevesQue et all (l2005i) 6 - IWalborn fc BladesI (Il997h . 7 - iHumphrevs et al.l (Il99'itl. 8 - IWaldron et al.l (|2004l '). 9 - 
iNieuwenhuiizen fc de Jagerl l|2000^ ■ 10 - Ide Jagej (|l998l ). 11 - ICrowther et all l|2006l ). 12 - iLobel et al.l (|2003l ) 
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Fig. 1. The behaviour of the factor F = 1 — L/L^^^ in the outermost meshpoint of the models along the sequence of 
homogeneous models. 
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Fig. 2. Evolutionary tracks of stars with Mo=60 and 120Mq in the Hertzsprung-Russel diagram for different assump- 
tions of the mass-los s prescription. For 60 Mq, the tracks for mass-loss prescription given by Eq. (4) with a=0.25 and 
mass-loss rates from Vink et al. ( 2000f) are shown. For 120 M©, tracks for Eq. (4) prescription with a=0.25 and 0.1 
are shown. The long-dashed broken line shows the Humphreys-Davidson limit. The star symbol marks the bending of 
the sequence of homogeneous models. 
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Fig. 3. Evolutionary tracks of stars in the Hertzsprung-Russel diagram. For 60 Af0, the track for a— 0.25 is shown; for 
120 Mq, tracks for a=0.25 and 0,5 are shown; for 200, 500 and 1001 A//©, tracks for a=0.25, 0,5, 0.8, 0.9 and 1.0 are 
shown. The dotted hne to the left shows the locus of homogeneous models. The dashed line indicates the Humphreys- 
Davidson limit (extrapolated for high Teg). The star symbol marks the position of a 133 M© star (indicating the 
bending point of the Z AMS). For o r ientat ion, the hatched rectangle shows the position of the Pistol star for the "low- 
luminosity" solution of iFiger et al. (1998) and the diamonds show the po sition of 77 Car, which is somewhat uncertain 
because of the ambiguity in the radius determination (jHillier et al.ll200ll ). 




t, yr 

Fig. 4. Stellar mass vs. evolutionary lifetime for a=1.0. Initial mass of stars is 1001, 843, 500, and 2OOM0. The dashed 
lines connect the loci of models with a hydrogen abundance in the convective cores of Xc =0.6 and 0.5. The star 
symbols at the curves indicate the central hydrogen exhaustion time (see also Table 1). 



1000 



\500 







_ 1 


1 1 


' 1 


1 1 1 


1 ' ' ' ' 1 - 




=0.6 






a=0.25 


Y=0.3 












0.4 


























^Y3=0.0 






I 1 1 1 1 I 1 1 1 1 I I 

10^ 2x10^ 3x10^ 

t, yr 



Fig. 5. Stellar mass vs. evolutionary lifetime for Q:==0.25-sequcnccs of stellar models. The dotted lines connect the loci 
of models with surface helium abundance 0.3, 0.4, 0.5, and 0.0. The dots at the curves indicate the time and mass 
when the hydrogen abundance in the convective cores Xc becomes 0.6(0.1)0. 
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Mass-loss rate vs. evolutionary lifetime for models with a— 0.25. The left border of the plots corresponds 
10'°^ yr, since mass-loss rates before this time are virtually the same as for the left extremes of the curves. 
The star symbols at the curves indicate the central hydrogen exhaustion time. Vertical bars show the ranges of the 



Fig. 6. 

to t = 



estimates of the mass-loss rate, obtained by means of q uantitative spectrosc opy for the Arches cluster (jNajarro et al 
2004 AN), Hcl emission stars in the Galactic center ([Martins et al.l 120071) . star- forming-regions in the Galaxy and 



the Magellanic clouds (AdK, A.de Koter, unpublished), and estimates for the Arches and Quintuplet stars obtained 
by radio-observations (AL and QL, respectively, .Lang et al. 2005.). F or the Quintuplet, only M for stars identified 
with stellar winds are shown. In the Arches cluster, Lang et ahl l 2005f ) identify stellar winds for all observed sources. 
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Fig. 7. The variation of L/L^dd in the outer meshpoint of the models along evolutionary sequences for stars with 
A/o=1001, 500, 200, 120, and 60 Mq, for a=0.25. The star symbols label models with hydrogen abundance in the 
convective core Xr — 0. 
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Fig. 8. The variation of the mass-loss rate versus Toff for the a— 0.25 sequence of stars (sohd curves). For comparison, 
the dotted curves represent estimates of the maximum M along the same tracks according to Eq. [SI 
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Fig. 9. The variation of the "wind performance number" along evolutionary tracks for a~0.25. 
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Fig. 10. The relation between stellar mass and mass-loss rate for the a— 0.25 case. Ticks along the tracks (in the 
order of decreasing mass) indicate models in which the surface He-abundance becomes = 0.4, the last model 
with Tee > 10000 K and Ys > 0.4, and the first model in which T^e > 10000 K again. The dots at the tracks 
mark the TAMS. The long-dashed line shows the mass-loss rate computed using the lVink et al.l (j200l[ ) mass- loss for- 
mula. Dotted, dot-dashed , and d ot-dot-dot-das hed li nes show mass-loss rates for hydro gen-deficient WR stars from 
de Ponder fc Vanbeveren ( 2003 ). Langer ([l989l), and Nelemans fc van den Heuvel ( 2001 ). respectively. Black squares 
represent estimates of stellar parameters of observed hydrogen-r ich WR stars based on model atmospheres that ac- 



count for iron- line blanketing and clumping ( Hamann et al. 2006[ ): open circles are estimates by the same authors for 



hydrogen-poor stars {Xh < 0.2) Asterisks are spectroscopic mass and M estimates for Galactic 0-stars, plusses are 
estimates for LMC 0-type stars, crosses are estimates for SMC stars (see Table [5] for references). The star symbol 
marks the position of the almost identical components of WR20a, currently the most massive star weighed in a binary 
system. Note its low M. Vertical bars have the same meaning as in Fig. [6l 
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Fig. 11. Evolutionary tracks of stars in the Hertzsprung-Russel diagram (solid lines, a=0.25 case). Dots along the tracks show positions of models with surface 
helium abundance Fs=0. 3(0. 1)0.9. Three dotted lines stretching from top to bottom connect ZAMS stars, the loci of models with surface He-abundance Yc = 0.4 
(models with higher Yc are conventionally considered as WR-stars), and TAMS stars. Broken long-dashed line shows the Humphreys-Davidson limit (linearly 
extrapolated at high L and Teff). The hatched rectangle shows the position of the Pistol star (see the text for comments). Also plotted are the positions of 
the most lum inous stars in the G alaxy, the LMC, and the SMC (Table [2]), the Hel-emission line stars in the Galactic center (jNajarro et al.lll997[ ). the Arches 



cluster stars ( Najarro et al. 2004( l. luminous WNLh stars in the Galaxy and the M a.gellanic clouds (AdK, A . de Koter, unpublished) and the most luminous 



stars in the R136 cluster, with spectral subtypes 03 If*, WN, 04If, and 03 III(f*) (jMassev fc Hunteij[l998l Table3). For AdK-stars, open circles are for the 
Galaxy, filled circles are for the SMC, crossed circles are for the LMC. See Appendix for the notes on particular stars. 
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Fig. 14. Initial-final mass relation for SMS and inferred outcomes of evolution: black hole formation for progenitors 
with mass M, ^ 250 Mq or Mj ^ 800 Mq, and pair- instability supernovae in the intermediate range of Mj. 



